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The interaction between an atom moving in a model double-well potential and the conduction 
electrons is treated using renormalization group methods in next-to-leading logarithmic order. A 
large number of excited states is taken into account and the Kondo temperature Tk is computed as 
a function of barrier parameters. We find that for special parameters Tk can be close to IK and it 
can be of the same order of magnitude as the renormalized splitting A. However, in the perturbative 
regime we always find that Tk J; A with a Tk ^ IK [Aleiner et ai, Phys. Rev. Lett. 86, 2629 
(2001)]. We also find that A remains unrenormalized at energies above the Debye frequency, o>Debye- 
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I. INTRODUCTION 

There is a great number of experimental anomalies ob- 
served by point contact [1-3], dephasing [4,5] and trans- 
port [6,7] measurements in metals at low temperature 
which have not been theoretically explained in a satisfac- 
tory way. Sensitivity on annealing [1,5], sample prepara- 
tion [5], electro- migration [2] support that these anoma- 
lies may originate from some kind of dynamical defects. 

The concept of two-level systems (TLS's), i.e. atoms 
or groups of atoms moving between two positions result- 
ing in two almost degenerate levels, has been introduced 
to explain the low temperature specific heat anomalies 
in metallic glasses [8]. In order to keep the level split- 
ting very low, tunneling between the two positions has 
been assumed. The TLS model has later been gener- 
alized by incorporating dissipation effects [9] and non- 
commutative couplings between the TLS and conduction 
electrons such as screening and assisted transition due to 
the metallic electrons [10,11], and a two-channel Kondo- 
like behavior (2CK) has been conjectured [11]. The es- 
timated Kondo temperature Tk, however, turned out to 
be too small. Taking a few of the higher levels of the 
atom into account increased Tjc considerably [12], how- 
ever, this approximation has been shown recently to be 
incorrect [13,14]: Most of the terms cancel out at energies 
above that of the few lowest levels if all excited states are 
included in the computation, and thus electron assisted 
tunneling results only in a negligible Tk for TLS's [13,14]. 
Similar results were obtained using the adiabatic approx- 
imation for an atom moving in a metallic host [15]. 

The physical reason of the above cancellation is that 
the tunneling takes place on a typical time scale r ~ 
l/u>D- Therefore electrons that are farther from the 
Fermi surface than ujd follow the motion of the defect in- 
stantaneously, and can be ignored when considering com- 
plicated correlations between the TLS and the conduc- 



tion electrons. In other words, the effective bandwidth 
of the conduction electrons is reduced from the Fermi 
energy Ep down to a value ~ u>d- 

In the present paper we study a similar model but 
with smaller or negligible barriers, far from the tunnel- 
ing regime, and keeping all excited states through our 
computations. 

We construct the perturbative scaling equations up to 
next to leading logarithmic order and show that the can- 
cellation found in Ref. [13] also extends to the splitting 
of the TLS, which remains unrenormalized down to the 
scale WDebye too. This supports again the picture of 'adi- 
abatic conduction band', proposed in Ref. [15]. 

In the most interesting cases the second level is just 
above or around the top of the barrier: In this regime we 
find a Kondo temperature in the range of Tk ~ 0.1 — IK 
using realistic parameters. For a special class of parame- 
ters the renormalized energy splitting between the lowest 
two levels is around Tk implying that such dynamical de- 
fects may give rise to some of the Kondo-like anomalies 
observed [1,3,6]. Our calculations also indicate that to 
obtain a measurable Tk one needs resonant scattering 
on the defect. 

In the present work we only study the perturbative 
region, and there we do not find convincing evidence of 
an observable two-channel Kondo behavior. However, 
one can show by performing analytical and numerical 
calculations using a more detailed model that there exists 
a regime outside the range of perturbative calculations 
where the two-channel Kondo behavior appears [16]. 

II. MODEL 

Our model consists of a particle (or collective coor- 
dinate z of the defect) with mass M (M ~ 50m p ; m p 
the proton mass) moving in some bare potential well, 
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Vbaie{z)- The potential Vbare(z) can be thought of as the 
potential resulting from the interaction with the neigh- 
boring ionic charges. However, as we shall see later, Vbarc 
is strongly renormalized due to the strong electron-defect 
interaction: In general, both Vbarc (z) and the electron- 
defect interaction are of the order of the Fermi energy 
(~ lOcV), but their sum, the effective potential is only 
of the order of the Debye frequency. Therefore Vbarc (z) 
has very little physical meaning. 

The non-interacting part of the Hamiltonian is given 
by: 
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The first term describes the motion of the particle in the 
potential V b ar e : b\ (n = 1,2,.. .; b\b n = 1 ) denotes 
the creation operator corresponding to a state of the par- 
ticle with energy e n and wave function ip n (z). The second 
term describes the conduction electrons, a pj(T being the 
creation operator of an electron with momentum p, spin 
tr, and energy e p . In the following we shall use a simple 
free electron approximation to describe the conduction 
band, and assume that the wave function corresponding 
to a p a is a simple plane wave ~ e lk r . The density of 
states for the electrons is g(e) « go{l + as/D ), where Qo 
is the value at the Fermi energy e = 0, Dq is a symmet- 
ric bandwidth cut-off of the order of the Fermi energy, 
and | a | < 1 accounts for possible electron-hole symmetry 
breaking due to band structure effects. 

We use a simple local interaction between the particle 
and the electrons: 

Hint = J dzdz p il>l(z)il) a (z)&(z p )TS>(z p ) 
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where ip a (z) is the electron field operator along the axis 
z of the defect motion (x = y = 0) and &(z p ) — 
bn tfnizp) is the particle field operator at the position 
of the heavy particle, z p . Uq describes a static scattering 
at the center while U describes the change in the scat- 
tering potential when the particle is displaced from the 
center. 

In the following we shall simply take Uq = 0. This 
simplification needs some explanation: In principle, the 
value of U should be precisely Uq in Eq. 2. One can, 
however, argue, that the electronic wave functions corre- 
sponding to the operators a p in Eq. 1 should be deter- 
mined when the particle is at the origin, and should there- 
fore already incorporate the effect of Uq, which could 
therefore be ignored. 

The previous argument is, however, not quite right. 
The consequences of the presence of Uq would be quite 
important. The term ~ Uq can be treated exactly for a 
single impurity [17]: Since Uq describes a static potential 



it can be incorporated in the electronic wave functions, 
which cannot be approximated by plane waves any more. 
Its most important effect is to renormalize the local den- 
sity of states qq in the s-channel as qq — > cos 2 (Jo) Qo, 
where So denotes the scattering phase shift induced by 
the scattering potential Uq = U . As we shall see below, 
experiments indicate that the value of U is rather large, 
and corresponds to almost resonant scattering. In other 
words, in this simple model with U = Uq at resonant 
scattering Sq « 7t/2, and the local density of states is 
strongly suppressed. As a result, the dimensionless cou- 
plings estimated below are reduced. This is a serious 
problem: One can show that for a potential scattering 
model with U = Uq this suppression is so strong that 
one is always in the weak coupling regime, and one can 
therefore never observe the two-channel Kondo behavior, 
while without the term Uq one can go beyond the Emery 
- Kivelson line and prove even analytically the existence 
of a regime where a two-channel Kondo behavior appears 
[16]. 

To resolve this problem, one has to go beyond approx- 
imating the particle by a simple potential scatterer and 
one must take into account the dynamics of the internal 
electronic degrees of freedom of the particle [16]. This 
analysis turns out to justify our approximation of setting 
U = in Eq. (2). 

To handle the interaction part of the Hamiltonian we 
follow Rcf. [11] and introduce spherical coordinates. As- 
suming that the motion of the defect is restricted in space 
(kpz p < 1, with kp the Fermi momentum), the domi- 
nant electron scattering occurs in the (I = 0, m = 0) and 
(I = l,m = 0) angular momentum channels, and we can 
neglect scattering channels with I > 1 [11]. H lrA then 
becomes: 
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where a\ l(J creates an electron with angular momentum 
I (m = 0), radial momentum k and spin a, and satisfies 
the anti-commutation relation {ffl^, dk'l'a'} = 27r 6(k — 
k') 8 aa > 5u> . We evaluated the dimensionless interaction 
matrix elements V^™ using spherical wave functions with 
momentum k « kp for the electrons and the exact wave 
functions ip n for the defect. The dynamics of the heavy 
particle is described in terms of the pseudo-fermions 6^ 
satisfying the constraint J2 n b\ib n = 1- 



III. RENORMALIZATION GROUP ANALYSIS 

To determine Tk we constructed the leading and next- 
to-leading logarithmic renormalization group (RG) equa- 
tions [12]. To this end one has to compute vertex and 
self-energy corrections to the pseudo-fermion propagator 
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Qnn'ij) = — (T T b n (T)b n i (0)) and the impurity-electron 
vertex function T™ , and then apply the relatively stan- 
dard machinery of multiplicative RG. The corresponding 
skeleton diagrams are shown in Fig. 2. [Note that the 
first order self-energy diagram (not shown in Fig. 2) only 
generates a local non-logarithmic and time-independent 
renormalization of the defect potential and can there- 
fore be entirely ignored.] In the RG procedure the band- 
width cut-off is reduced D — > D + dD and the couplings 
are simultaneously changed to keep physical quantities 
invariant. In the present case a matrix version of the 
multiplicative RG must be employed [18]. 

In this work we take also into account the energy de- 
pendence of the local density of states. This may also de- 
pend on the angular momentum channel /, Qi(e)/gi(0) = 
1 + ceie/Do, where the slope a; and value gi(0) of the 
density of states is usually different for different values of 
I. The energy e is measured from the Fermi energy and 
only the linear term has been kept in gi. 

Let us first focus on the case of gi = o — gi=i — go and 
a n = ai = a. Then the second order self-energy cor- 
rection of the pseudo-fermion propagator, e.g., contains 
both logarithmic and non-logarithmic corrections, and is 
given by the following expression: 
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where we introduced a matrix notation, V™ — ► V_ nn , 
tr{...} denotes the trace in the electronic angular mo- 
mentum indices, and the constants f(a) and g(a) are 
given by 



/(a) = -a 2 (ln2-l) + 21n2, 
g{a) = 1 -ln2 + a 2 (i - In 2) . 
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The factor two in Eq. (4) is due to the spin degeneracy 
of the electrons. The effect of the non-logarithmic terms 
proportional to / and g is to strongly renormalize the 
eigenstates and eigenfunctions of the double well poten- 
tial by replacing the the heavy particle Hamiltonian by 
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where the matrix A nm is given by 

A nm = e n S nm - 2 tr{X n ":r m } [Do /(a) - e n g(a)] . 

h 

(8) 

The effect of the static (w-independent) non-logarithmic 
terms in Eqs. (4) can be taken into account by diagonal- 
izing A nm through a unitary transformation 



(UAU*) nm = l m S nm , 

ynm _^ (Xjyjjjyim ^ 
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where the e„'s denote the renormalized values of the 
heavy particle energies. The effective Hamiltonian corre- 
sponding to A nm generally contains non-local terms too, 
but the largest terms actually turn out to be simple local 
corrections to Vbare- Therefore, in the rest of the paper we 
shall drop these non-logarithmic self-energy corrections, 
and assume instead that we can model the entire effec- 
tive Hamiltonian by a simple local quasi one-dimensional 
symmetrical square potential with a barrier in the middle 
and infinite walls (see Fig. 1). We determine the corre- 
sponding wave functions by solving simple transcendental 
equations and then use these wave functions to compute 
the interaction matrix elements y_ nm . 

Note that it is only the renormalized Hamiltonian that 
can be measured and has therefore physical meaning. 
Note also, that both Vbaro and the the corrections are 
usually of the order of the Fermi energy, however, their 
sum is usually of the order of the Debye energy only. 



SB 


z-- 

-' 


=0 

I / 






/ > 



2a^a 



B 



FIG. 1. The simple square-shaped double-well potential 
used: Vb is the height of the barrier, a ~ 0.1 A is the width of 
the wells and 2aj,a ~ 0.4 — 0.5A is the barrier width. Inset: 
Effective scattering potential for the electrons. The shaded 
potential corresponds to the defect sitting at the center. The 
white potential is due to the displacement of the defect. 

The interaction vertex also contains a non-logarithmic 
part. Up to second order in V_ the dimensionless vertex 
function can be calculated as 

7 nm (a>) = V nm - q ^ {ynkykm + yfrmynk^ 
e k <D 
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The effect of the non- logarithmic terms in (11) can be 
taken into account by renormalizing the bare vertex as 
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(12) 
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This transformation must be constructed order by or- 
der, and it sums up systematically all higher order non- 
logarithmic vertex contributions. 

After performing the transformations above the ex- 
pressions of Q" 1 and 7 simplify considerably and contain 
only logarithmic terms in Dq: 

[S^WM = w - S nm € n (13) 

■nk -km , , _ , Do 
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+ 2 2^ tr{j/ F }(w - e fc ) ln 



7 „ m(w) = £»m _ y ln _D0 ^nft^fimj (14) 



The remaining logarithmic terms in Eqs. (13) and (14) 
can be summed up using a generalized version of the 
multiplicative RG [18], leading to the following RG equa- 
tions: 
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Here x = ln(_D /-D) is the scaling variable, and D de- 
notes the initial value of the cut-off. As the Hamiltonian 
is diagonal in the spin, each closed electron loop results 
in a multiplicative factor of 2 corresponding to A// = 2 
conduction electron channels. 

The renormalization group equations for the energies 
e n are somewhat more complicated. In particular, the 
RG generates off-diagonal terms to the heavy particle 
Hamiltonian: 
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where e„,e m < D, and dx — \n[D/D'}. Therefore, in 
each RG step we diagonalize the self-energy by construct- 
ing a new set of defect states with renormalized energy 
eigenvalues, and express all couplings in this new basis: 
(p n = ^2 m U mn ip m . Note that Eqs. (15) are invariant 
under this unitary transformation. 

When the reduced cutoff reaches the renormalized en- 
ergy eigenvalue of some defect level, the dynamics of that 
level is frozen out, and we therefore drop it in the follow- 
ing RG steps. Usually only a few levels remain active 
in the region where the relatively weak initial couplings 
become comparable to the stronger ones. In many cases 
only one level remains by this time, meaning that no 
Kondo effect occurs at all. 




FIG. 2. Diagrams used to construct the renormalization 
group equations: Dashed lines correspond to the defect atom, 
solid lines to the electrons. 

In the original tunneling model [11] only the states 
n = 1 , 2 were kept giving rise to two initial dimension- 
less coupling constants, v x and v z , corresponding to the 
diagonal and off-diagonal terms in indices n and n'\ 



v* = \{V£-V£-V™ + V??). 



(17) 
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We identify Kondo temperature Tk with the energy scale 
at which the initially small dimensionless coupling v x 
reaches about the half of its fixed point value, v x ~ 0.1 
[19], and we keep the last two levels active even if their 
separation is larger than the running cut-off D. Be- 
low the Kondo temperature the perturbative RG breaks 
down and more efficient methods are needed. If only 
two levels are kept then Bethe-ansatz [20] results for the 
two-channel Kondo problem could be applied. 

The bare values of v z can be estimated from ultrasound 
data [11] and for alloys with stronger couplings they cor- 
respond to v z ~ 0.2 and thus Uqo ~ 3 for typical param- 
eters. This rather large value of Uqo implies resonant 
scattering on the atom [21]. Nonetheless, even Ugo <~ 3 
turns out to be small enough so that the couplings V nn 
corresponding to it are still in the perturbative regime: 
Restricting our discussion to n = 1, 2, the largest matrix 
elements are those where both the electron and defect 
parities are changed and are about v z ~ 0.2. The re- 
duction from Uqo ~ 3 is due to a factor v z ~ kpd U go, 
with d = 2a(l + a^) ~ 0.4 A the width of the potential 
well, and hp 1 ~ lA the Fermi momentum. Couplings 
with conserved parities such as v x are even smaller, and 
are further reduced by the Gamow factor in the case of 
a large barrier in the tunneling regime [11]. 



IV. RESULTS 

In this work we used the following procedure: We first 
diagonalized the effective heavy particle Hamiltonian to 
obtain the renormalized defect energies e n and the renor- 
malized couplings V numerically. Then we summed up 
logarithmic terms by performing the RG steps described 
in the previous section. 
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Firstly, to determine the fixed point structure of the 
flow equations and to determine the Kondo temperature 
we continued the RG even after reaching the first excited 
state (where the dynamics of the defect must be entirely 
frozen due to the energy splitting) and verified that then 
the couplings indeed flow to the two-channel Kondo fixed 
point as expected. 

The leading logarithmic scaling equations for many 
levels were investigated by Alcincr et al. in the tunneling 
regime [13]: They correspond to the first two diagrams in 
Fig. 2 and to the second order terms in Eq. (15). Aleiner 
et al. have shown that in leading logarithmic order the 
logarithmic contribution of the excited states cancels out 
due to a sum rule, related to the approximate complete- 
ness of the defect wave functions. This reduces the ef- 
fective cutoff to the order of the Debye temperature and 
thus the essential renormalization of the couplings occurs 
in the region where only few excited states remain active. 

We find that this is also true for the renormalization 
of the tunneling rate. Fig. 3 shows the renormalization of 
the couplings, v x , and the tunneling rate, A12 = §2 — ?i, 
computed by solving leading and next-to-leading loga- 
rithmic scaling equations numerically. Though we could 
not find a simple sum rule similar to that of Ref. [13], 
we found that both v x and A12 remain unrenormalized 
throughout a large energy region, down to an energy scale 
of the order of the Debye temperature ~ cjjj. This is 
in agreement with the 'adiabatic conduction band' pic- 
ture of Kagan and Prokof 'ev, [15] who argued that, since 
the tunneling takes place during a time scale ~ I/Vd, 
electrons with energies larger than u>u follow the defect 
instantaneously and therefore do not contribute to the 
logarithmic singularities summed up by the renormaliza- 
tion group. Therefore the drastic changes in the effective 
couplings appear only when a few levels remain. 

For typical bare values of the splitting A = A21 = 
E2 — E% > IK, the renormalization of A turned out to 
be less than about 25%. This must be opposed to the 
case where only two levels are kept from the beginning, 
and therefore a strong power-law renormalization occurs 
even at high energies [11]. Our results agree with those 
obtained by means of an adiabatic approximation where 
the renormalization occurs only when the cutoff is al- 
ready in the range of the few lowest energy levels [15]. 
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FIG. 3. The RG flow of A 12 and v x in next-to-leading log- 
arithmic order for a = 0.4, d = 0.5A and Vb = 100K. Only 
curves corresponding to goU < 2.0 are shown, since for larger 
Qoll the weak coupling approximation is not valid. 

It is very instructive, how Tk depends on the number 
of defect states (see Fig. 4) . Such a calculation has been 
first performed by Aleiner et al. [13]. The first few states 
result in a strong increase in Tk, but then Tk suddenly 
drops, and the states n > 6 practically don't produce 
any further changes. In the following only the 'satura- 
tion values' of Tk will be reported that we obtained by 
eventually continuing the RG flows below the first excited 
state if A 2 i > Tk- 
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FIG. 4. The RG flow of the couplings in leading logarithmic 
order and in next-to-leading logarithmic order. Inset: Kondo 
temperature as a function of the number of defect states, for 
Uqq = 2.1, d = 0.5A and Vb = 100K . 



It is very important to determine the range of validity 
of the perturbative scaling analysis we use. Fortunately, 
there is a few non-perturbative results available for the 
TLS model that can be used to achieve this goal. In the 
limit of small A12 and v x the scaling equations take a 
simple form [22]: 
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where the scaling exponents yt and y x of the dimcnsion- 
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These approximate exponents are compared to the exact 
ones in Fig. 5: The exponent y x of the assisted tunneling 
is surprisingly well approximated by y x G in the range 
v z < 0.25, and therefore we expect that the estimate 
of the Kondo temperature is reliable in the range v z < 
0.25 — 0.3. However, yf- G underestimates the value of 
yf- G and therefore overestimates the renormalization of 
A12, for which the RG results should be trusted only in 
the range v z < 0.22. 
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FIG. 5. Comparison of the exact scaling exponents yt and 
y x with their approximate values y^ G and y^ G . 

Fig. 6 shows the barrier height dependence of the 
Kondo temperature: Tk increases with increasing bar- 
rier height up to a certain value of Vg and then it drops 
suddenly, while the splitting A21 decreases continuously. 
Initially, increasing the barrier height concentrates the 
wave functions of the first two states more and more 
around the potential minima and this results in an in- 
crease of the coupling v z , and thus a gradual increase in 
Tk- This tendency is, however, suddenly reversed once 
the barrier reaches the first level, and one enters the tun- 
neling regime (Vb > 150K for the parameters of Fig. 6): 



There the assisted tunneling v x , responsible for the gen- 
eration of the Kondo effect, decreases exponentially, and 
both A12 and T% decrease dramatically. 
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FIG. 6. The 
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The Kondo temperature as a function of the height 
of the central barrier, for a = 0.4 and M — 50m p . The 
renormalized level splittings are plotted as well. The shaded 
part of the figure indicates the region where the estimate of 
the renormalized A12 is unreliable. 



Tk is also very sensitive to the width of the entire 
potential well d. In Fig. 7 we show the G?-dependence of 
Tk for a defect that has a relatively low barrier height 
and is not in the tunneling regime. As d is decreased from 
d = 0.5A, the energy levels are shifted to higher values 
and also the initial values of the coupling v z decrease 
while v x does not change dramatically. Thus the width of 
the well is a fundamental parameter and the Kondo effect 
can only occur if the room for the particle considered is 
anomalously large in one direction. For defects in the 
tunneling range Tk decreases with increasing d since the 
bare value of v z only slightly increases while the assisted 
tunneling v x is exponentially suppressed in the tunneling 
regime. 
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FIG. 7. The Kondo temperature and the renormalized level 
splittings as a function of the total width d of the potential 
well for q = 0.4 and M — 50m p . The shaded range of the 
figure indicates the region where the estimate of the renor- 
malized A12 is unreliable. 

In Fig. 8 we show the dependence on the strength of 
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the defect-electron interaction U. Again, Tk increases 
continuously with increasing U, since all couplings gener- 
ating the Kondo effect increase. On the other hand, The 
renormalized A decreases continuously since the scaling 
exponent yf" also decreases with increasing U. 
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FIG. 8. The Kondo temperature and renormalized level 
splittings as a function of the scattering strength U. The 
shaded part of the figure indicates the region where the RG 
estimate of the renormalized A12 is unreliable. 

Finally, let us discuss the electron-hole asymmetry 
dependence of T#. For Uq = 0, qi = q = Qi = i and 
a = a a = ai we find no significant a-dependence of 
the renormalized splitting and Tk- In general, however, 
the local density of states is different for the / = and 
I = 1 channels and thus both qq and a depend on I. 
A non-zero potential scattering term Uq ^ in Eq. (2), 
e.g., would have a twofold effect: (1) It would change the 
local density of states in the / = scattering channel, 
and (2) would generate a different electron-hole symme- 
try breaking for the I = and I = 1 channels [17]. [Note 
that in Ref. [23] Uq has not been set to zero and has 
been treated inappropriately: Instead of incorporating it 
in the parameters Qq and «o we included it in the RG 
equations as a coupling constant. Though the results ob- 
tained are qualitatively similar to the ones we obtained 
below, this procedure is wrong, since mass terms like Uq 
have to be treated always differently from the coupling 
constants in the multiplicative RG [17].] 

It is easy to generalize our previous calculations to the 
case cto 7^ a\ and qi=q ^ Qi=i- The /-dependent den- 
sities of state can be treated by defining the fermion 
fields slightly differently. However, for o.q ^ oc\ the 
electron-hole symmetry breaking generates strong non- 
logarithmic corrections to the bare coupling constants. 
In Fig. 9 we plotted the dependence of Tk on ao while 
keeping a± zero and qi—q — Qi—\. We find that for typi- 
cal parameters we used earlier Tk can change about an 
order of magnitude. This renormalization can have even 
more dramatic effects for very small Tr-'s where Tk can 
change several orders of magnitude due to changes in the 
value of the electron- hole asymmetry (see inset of Fig. 9). 
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FIG. 9. Electron-hole asymmetry dependence of Tk and 
the splitting A12. Inset: the effect is more dramatic for very 
small Tk 's. 

Let us conclude this Section with a general observa- 
tion: 

Increasing the mass (the M/m p ratio), the energies of the 
levels are scaled down. If the central barrier is similarly 
reduced then Tk becomes smaller by the same ratio and 
thus A/Tk is not affected, as the dependence on the high 
energy cutoff D is weak. 



V. CONCLUSIONS 

In this paper we performed a thorough scaling analy- 
sis of a dynamical defect coupled to the conduction elec- 
trons, by taking into account all the excited states of the 
defect, and constructing the next to leading logarithmic 
scaling equations. We also discussed how to take into ac- 
count large non-logarithmic terms that renormalize the 
bare parameters of the theory. In the perturbative regime 
our calculations confirmed the results of Ref. [13]: The 
electron-defect couplings remain unrenormalized down to 
an energy scale of the order of the Debye temperature 
D ~ WDebye due to a cancellation pointed out in Ref. [13]. 

In the leading logarithmic approach of Ref. [13] the en- 
ergy levels of the impurity are unrenormalized, and the 
above next to leading logarithmic analysis is needed to 
compute the renormalization of the defect states. We 
found that, similar to the defect-electron coupling, the 
splitting of the first two defect states remains unrenor- 
malized down to WDcbyc- In contrast to the expectations 
of Ref. [24] , where the generation of a large splitting has 
been predicted, we found that the renormalization of the 
defect states remains small if one takes into account all 
defect states. 

Our results perfectly agree with the adiabatic picture 
of Kagan and Prokof'ev [15]: Electrons with excitation 
energies |e| > WDebye follow the defect instantaneously, 
and therefore the role of the excited states is only reduc- 
ing the effective electronic cut-off from the order of the 
Fermi energy down to the energy scale WDebye- 
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We also analyzed the range of validity of our approach 
by using some analytical results in the strong coupling 
regime. We found that the perturbative RG predicts the 
Kondo temperature correctly for surprisingly large cou- 
plings, however, it tends to overestimate the renormal- 
ization of the splitting A. 

We solved the RG equations for initial couplings ob- 
tained by changing various parameters of the double 
well potential. We always found a region where Tk 
and the renormalizcd A12 were comparable and there- 
fore one should be able to observe the Kondo anomaly 
in various measurements. Our calculations show that to 
have Tk ~ A12 in the measurable range one needs a 
defect that (i) has anomalously large room to move in 
one direction, (ii) is close to but not yet in the tunnel- 
ing regime (Hi), which has a large scattering strength, 
Uga > 2.5, implying resonant scattering on the parti- 
cle. The best candidates are therefore atoms with res- 
onant d or f-scattering at the Fermi energy or possibly 
small groups of atoms [25] . Thus dynamical local defects 
with special realistic model parameters can explain the 
Kondo-like anomalies observed in some experiments for 
T > Tk [1,3,6]. Note that in amorphous materials the 
positions of the defect /- or ci-orbitals have usually a dis- 
tribution, and therefore many of the defects may have 
resonant scattering at the Fermi energy. 

However, to explain the two-channel Kondo scaling re- 
ported in Rcf. [26], one would need a renormalizcd A12 
that is less than Tk at the energy scale D ~ Tk ■ Unfortu- 
nately, our calculations are only of logarithmic accuracy, 
and therefore we cannot decide if the ratio A 12 /Tk is 
small enough to display a clear two-channel Kondo scal- 
ing. In the perturbative regime (i.e. the regime where 
where our perturbative RG works), our results seem to 
indicate that Tk can be comparable to, but is somewhat 
smaller than A12, and is thus too small to result in an 
observable two-channel Kondo behavior. However, one 
can prove both analytically and with NRG calculations 
the existence of a regime outside the range of perturba- 
tive RG [16] where Tk can be larger than A12. This is 
also indicated by our perturbative results. 

We have to emphasize that for a defect with resonant 
scattering in a disordered environment many of our ap- 
proximations (free electron model, simple delta scatter- 
ing on the defect, usage of a local effective potential 
etc.) are questionable, and the estimates of the vari- 
ous couplings are therefore inadequate. Also, in reality, 
the dynamical defect is not formed by a single atom, 
but it is probably rather a group of atoms (disloca- 
tion, e.g.) that is responsible for the defect dynamics. 
Though our general conclusions probably do not depend 
on these approximations, Tk, e.g., depends exponentially 
on the coupling constants and more sophisticated mod- 
els/calculations would be needed to give a quantitative 
estimate of Tk and the ratio A^/Tk- 

It remains an open question whether for larger cou- 



plings (i.e., defects with resonant scattering at the Fermi 
energy) a larger Tk ~ 5 — 10K can be reached, as required 
to explain the experiments for Rcfs. [1-3]. 
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